EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



X 



CERN-PH-EP/2004-010 
18. February 2004 



Measurement of the Strange Spectral 
Function in Hadronic r Decays 

o 

cN ■ The OPAL Collaboration 

(N 

Abstract 
> 

^^ . Tau lepton decays with open strangeness in the final state are measured with the OPAL detector 

^D I at LEP to determine the strange hadronic spectral function of the r lepton. The decays r~— > 

(Kvr)^!/^, (KTnT)~i'r and (KTnnT)~ Vr with final states consisting of neutral and charged kaons and 

'!::j- ■ pions have been studied. The invariant mass distributions of 93.4% of these final states have been 

^^ , experimentally determined. Monte Carlo simulations have been used for the remaining 6.6% and 

>^ I for the strange final states including rj mesons. The reconstructed strange final states, corrected 

for resolution effects and detection efficiencies, yield the strange spectral function of the r lepton. 

C^. The moments of the spectral function and the ratio of strange to non-strange moments, which are 

(— I I important input parameters for theoretical analyses, are determined. Furthermore, the branching 

fractions 5(r" -^ K'TT^Ur) = (0.471 ± 0.064stat ± 0.022sys)% and B{t- -^ K-vr+vr-z^^) = (0.415 ± 

0.059stat i O.OSlgys) % have been measured. 
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1 Introduction 



The r lepton is the only lepton heavy enough to decay into hadrons. A comparison of the inclu- 
sive hadronic decay rate of the r lepton with QCD predictions allows the measurement of some 
fundamental parameters of the theory. The inputs to these studies are the spectral functions that 
measure the transition probability to create hadronic systems of invariant mass squared s = m?. 

The energy regime accessible in r lepton decays can be divided into two different regions: the 
low energy regime which has a rich resonance structure where non-perturbative QCD dominates; 
and the high energy regime near the kinematic limit, s = m^ = (1.777 GeV)'^, where perturbative 
QCD dominates. The high energy regime in r lepton decays provides an environment where the 
strong coupling constant ag can be measured [1-6] because the perturbative expansion converges 
well and non-perturbative effects are small. The measurement of the non-strange spectral function 
of hadronic r lepton decays [7-9] has provided one of the most accurate measurements of as, and 
some very stringent tests of perturbative QCD at relatively low mass scales [10]. The spectral 
function of strange decays allows additional and independent tests of QCD and a measurement of 
the mass of the strange quark [1, 11-13]. 

The strange spectral function of the r lepton is obtained from the invariant mass spectra of hadronic 
r decay modes with net strangeness. The relevant decay channels are listed in Tabled The only 
contribution to the one meson final state^ r~ -^ K'Ur is not reconstructed in this analysis. (Kvr)" 
and (K?7)~ contribute to the two-meson final states. They have a total branching fraction of 
1.367%. The decay r~ -^ K~7r°i/^ with a branching fraction of 0.450% can be fully reconstructed. 
The T~ -^ K^TT'Vr channel {B = 0.890%) consists of two decay chains: t~ -^ Kf^ir'Vr (50%) and 
r~ -^ Kg7r~z^T- (50%) where the Kg decays to two charged (neutral) pions in roughly 68% (32%) 
of all cases. Only those final states where the Kg decays to two charged pions are considered. The 
same applies to r" -^ K^Tr'TT^Ur {B = 0.370%). This and r" -^ K-tt+tt-z/^ {B = 0.280%) together 
form the two most important contributions to the three-meson final states. Finally, the four-meson 



^In order to simplify the text we refer only to the decays of the negatively charged r lepton. Simultaneous 
treatment of the charge conjugate decay is always implied. 
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Table 1: Overview of all channels with net strangeness and their branching fractions [14]- 



final state (KTrvrTr)^ is detected via tlie decay r^ -^ ¥>r ir^ tt^ -k^ v^ ■ In addition, the decay K~r]7r^i'T 
contributes with 0.029%. Hence, 93.4% of ah decay channels of the multi-meson final states with 
open strangeness, (Kvr)^, (Kvrvr)^ and (KvrTTTr)", were measured. This paper describes the selection 
of these dominant channels and the measurement of their invariant mass spectra using data collected 
with the OPAL detector during the LEP-I period from 1991 to 1995. The remaining 6.6% and the 
final states including rj mesons are taken from Monte Carlo simulation. The spectral function is 
then determined from these spectra and the spectral moments are calculated. 

From an experimental point of view, one of the key issues of this analysis is the separation of 
charged kaons and pions via the measurement of energy loss in the OPAL jet chamber in the 
dense environment of multiprong r lepton decays. Substantial improvements have been achieved 
compared to previous publications [15]. In particular, these improvements have made it possible 
to obtain a reliable dE/dx measurement in an environment where three tracks are very close to 
each other. The reconstruction of neutral pions is based on the study of shower profiles in the 
electromagnetic calorimeter. Furthermore the identification and reconstruction of r lepton decays 
with Kg — > TT~^TT~ has been achieved with high efficiency and good mass resolution. 

The outline of this paper is as follows: Section [2 gives a short description of the OPAL detector 
concentrating on those components which are important for this analysis. In addition, the r lepton 
selection and the Monte Carlo samples used are discussed. Section |S1 continues with a discussion 
of the experimental aspects of this work. The selection of the strange hadronic r lepton decays is 
described in Section |11 In Section IHl the results for the branching fractions, the strange spectral 
function and the spectral moments are presented together with a discussion of the systematic 
uncertainties. The results are summarized in Sectional 



2 Detector and Data Samples 

2.1 The OPAL Detector 

A detailed description of the OPAL detector can be found elsewhere [17]. A short overview is 
given here of those components that are vital for this analysis. Charged particles are tracked in 
the central detector, which is enclosed by a solenoidal magnet, providing an axial magnetic field 
of 0.435 T. A high-precision silicon micro-vertex detector surrounds the beam pipe. It covers the 
angular region of | cos 9\ < 0.8 and provides tracking information in the r — cp direction^ (and z from 
1993) [18]. The silicon detector is surrounded by three drift chambers: a high-resolution vertex 
detector, a large-volume jet chamber and z-chambers. 

The jet chamber measures the momentum and energy loss of charged particles over 98 % of the 
solid angle. It is subdivided into 24 sectors in r — (j), each containing a radial plane with 159 anode 
sense wires parallel to the beam pipe. Cathode wire planes form the boundaries between adjacent 
sectors. The 3D-coordinates of points along the trajectory of a track are determined from the sense 
wire position, the drift time (r — (p) and a charge division measurement (z) on the sense wire. The 
combined momentum resolution of the OPAL tracking system is cr-p/p'^ ~ 1.5 • 10^'^ GeV^ . From 
the total charge on each anode wire, the energy loss dE/dx is calculated and used for particle 
identification. This measurement provides a separation between pions and kaons of at least 2a in 
the momentum range relevant for this analysis (3 GeV < p < 35 GeV) . 

Outside the solenoid are scintillation counters which measure the time-of-flight from the interaction 



^In the OPAL coordinate system the x-axis points to the center of the LEP ring. The z-axis is in the e beam 
direction. The angle 9 is defined relative to the z-axis and cj) is the azimuthal angle with respect to the x-axis. 



region and aid in the rejection of cosmic events. Next is the electromagnetic calorimeter (ECAL), 
which, in the barrel section, is composed of 9440 lead-glass blocks, approximately pointing to 
the interaction region, and covering the range |cos^| < 0.82. Each block has a (10 x 10) cm^ 
profile with a depth of 24.6 radiation lengths. The resolution of the ECAL in the barrel region, 
including the effects of the approximately one radiation length of material in front, is ue/E = 
^(0.16)2GeV/^ + (0.015)2. 

The hadron calorimeter (HCAL) is beyond the electromagnetic calorimeter and is instrumented 
with layers of limited streamer tubes in the iron of the solenoid magnet return yoke. The outside 
of the hadron calorimeter is surrounded by the muon chamber system, which is composed of four 
layers of drift chambers in the barrel region. 



2.2 Selection of r Lepton Candidates 

For the selection of r lepton candidates, the standard r selection procedure described in [16] 
is used. The decay of the TP produces a pair of back-to-back highly relativistic r leptons. Their 
decay products are strongly collimated and well contained within cones of half-angle 35°. Therefore, 
each r decay is treated separately. In order to have a precise and reliable dE/dx measurement 
and to avoid regions of non-uniform calorimeter response, this analysis is restricted to the region 
|cos^| < 0.68. To reject background from hadronic events, exactly two cones are required and 
a maximum of six good'^ tracks in the event is allowed. This background is further reduced by 
requiring that the sum of the charges of all tracks in each individual cone is ±1 and the net charge 
of the whole event is zero. A total of 162 477 r cone candidates survive these selection criteria with 
an estimated non-r background fraction of 1.5 %. 



2.3 Simulation of Events 

The r Monte Carlo samples used consist of 200 000 r pair events generated at \/s = rago using 
KORALZ 4.02 [19] and a modified version of TAUOLA 2.4 [20]. Modifications were necessary 
because all four- and five-meson final states with kaons (signal as well as background) are missing 
in the standard version. Since the resonance structure of these channels is poorly known, only 
phase space distributions of these final states were generated. In addition, the resonance structure 
of various final states was modified to give a better description of the data [21,22]. The branching 
fractions of the decay channels with kaons are enhanced in this sample so that it comprises roughly 
a factor of ten more r decays with kaons than expected from data. The Monte Carlo events are 
then reweighted to the latest branching fractions given in [14], which are used throughout the 
selection procedure. The Monte Carlo events were processed through the GEANT OPAL detector 
simulation [23]. 

The non-T background was simulated using Monte Carlo samples that consist of 4 000 000 qq events 
generated with JETSET [24], 574000 Bhabha events generated with BHWIDE [25], 792 000 n- 
pair events generated with KORALZ [19] and 1755 000 two-photon events using PHOJET [26], 
F2GEN [27] and VERMASEREN [28,29]. 

^A good track has a minimum number of 20 hits in the jet chamber, a maximum |d()| of 2 cm, a maximum \za\ 
of 75 cm, at least 100 MeV transverse momentum and a maximum radius of the first measured point on the track of 
75 cm. 



3 Identification of Hadrons in r Final States 



3.1 Energy Loss Measurement in r Decays with Three or More Tracks 

A crucial part of this analysis is the identification of charged kaons via energy loss measurement in 
the jet chamber. Since this is the only means of distinguishing between charged pions and kaons in 
OPAL, a very good understanding of the effects present in the multi-track environment in r lepton 
decays is vital for any analysis that requires particle identification. 

The high Lorentz boost (7 ~ 25) of the r lepton results in its decay products being contained in 
a narrow cone with a typical opening angle of 5°. In those cases where the final state consists of 
more than one track, the dE/dx measurement is known to be no longer reliable [15]. A systematic 
shift in the dE/dx distribution is observed which leads to a misidentification of charged pions as 
kaons and thus to a reduced sensitivity in those cases where particle identification is required. The 
reason is explained in the following text. 

When a charged particle passes through the jet chamber, it produces an ionization cloud. This 
cloud travels at a constant speed of vd ~ 53iim/ns in the homogeneous drift field to the sense 
wires, where it produces a signal pulse. The pulse is integrated over a time At = 200 ns, which 
corresponds to a drift distance of ~ 1 cm. The integral is proportional to the energy loss of the 
particle in this particular drift volume. The hits used for the energy loss measurement are the same 
as used in the track reconstruction, but they have to fulfill additional quality criteria [30]. In the 
following, they are called ^ dE/dx hits'. 

If an additional track passes through the same sector, a second pulse will be created which overlaps 
with the tail of the first one (Figure^. The contribution of the tail is determined by extrapolation 
of a reference pulse [30] which is normalized to the integral over the first pulse. The pulse used in 
all previous OPAL publications overestimates this tail contribution. The observed deviation in the 
energy loss measurement as a function of the distance between these two pulses (At) is shown in 
Figure 121 This plot is obtained using tracks in r data that have 'first hits' and 'second hits', and 
which have a momentum greater than 3GeV. Such tracks arise when a second track traverses the 
same sector of the drift chamber and the two tracks have at least one intersection point in r — (/> 
within the sensitive volume. Thus, for a given track, some of the hits do not need tail subtraction 
('first hits') and some hits do need tail subtraction ('second hits'). Using 'first hits' only (given 
there are at least 20 of them) (d£'/d2;)exp is calculated. The 'second hits' of the same track were 
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Figure 1: Illustration of the tail subtraction procedure, (a) shows one half-sector of the OPAL 
central drift chamber with three tracks. The arrow indicates the drift direction of the 
ionization cloud following a path given by a Lorentz angle of 20°; (b) illustrates the signal 
seen on the sense wire for two successive pulses. 
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Figure 2: Measured energy loss normalized to the expectation as a function of time difference be- 
tween two measured hits (^200 ns = 1.054 cnij. The observed deviations using the standard 
correction are shown for (a) 'second hits' and (b) 'third hits'. The error band reflects the 
la error band of the parametrization. (c) shows the same distribution after all corrections 
have been applied. For better visibility, the distribution for 'third hits' is shifted by —0.1 
in the plot. 



then used to obtain (di?/dx)nieas which is analyzed as a function of the drift time difference Ai 
between two successive hits. The left-right ambiguity, which in principle exists if two tracks pass 
on either side of the anode plane, is treated properly here, since the correction is applied based on 
measured drift times and not on the spatial separation of the tracks. 

For two hits as close as At = 200 ns, the observed deviation is of the order of 10% of the measured 
dE/dx. In the region between 400 ns < At < 900 ns the measured dE/dx is slightly overestimated. 
The standard correction is finally switched off for pulses with a drift time difference of more than 
1000 ns which produces the structure shown in Figure Efa). The deviation as observed for 'third 
hits' can be seen in Figure |21^b). When weighted with the At distribution of all measured hits this 
gives on average the correct dE/dx for high multiplicity events. In multiprong r lepton decays 
this At distribution differs from the average in the sense that the distribution peaks at smaller 
At. Therefore, the energy loss measurement tends to be smaller than it should be. To avoid this 



problem, previous analyses [15] exploited the dE/dx information only of those tracks closest to 
the anode plane to classify the r decay mode. This however significantly reduces the number of 
identified decays. 

For this analysis, a new reference pulse has been developed that avoids the shortcomings of the 
standard one. In addition, a parametrized pulse shape was used instead of a binned one to avoid 
artifacts like the dip at At ~ 500 ns. The new reference pulse is of the form 

/ , .At, ,, ,n , (At)2 

-Pnorm= UiAtexp( ) + paiAtf exp{- 



P2 Pi 

/a /At, ,^ ,2 / (At)^ 

+ U5Atexp( +p7(At^exp-^ — '- 

\ P6 P8 

+ pl+P2At (1) 

with two terms to describe the short-range and the long-range part respectively plus a linear 
contribution. The pi are parameters that are optimized for the multi-track environment in r lepton 
decays. The correction is applied to all hits. The effect of the correction described above can be seen 
in Figure [21^c) for 'second hits' and 'third hits'. Apart from this normalization correction, a further 
bias reduction is obtained by also correcting the shape of the reference pulse depending on the charge 
deposited by the preceding pulse. In general, the new reference pulse shows a steeper rise at low At 
and a lower tail to avoid overestimation of the tail subtraction for subsequent pulses as explained 
above. As a result of this procedure, a dE/dx bias reduction to ±1% in (d-E/dx)meas/(d-E/dx)exp 
has been achieved. As the procedure is applied iteratively, i.e. the first pulse is used to correct the 
second, the first and corrected second pulse are used to correct a possible third pulse and so on, 
the chosen reference pulse is valid for any jet topology [30]. 

If a track is close to the anode or cathode plane in the jet chamber, the drift field is no longer 
homogeneous and one observes a deviation in the measured dE/dx of the order of 3%. Corrections 
for this effect are determined using Z^ -^ ^~M^ events. 

When a signal is measured at a sense wire, an induced signal at the neighboring wires is also present. 
This effect depends on the momentum (or curvature) of the track and on cos 9. Corrections have 
been determined using |U-pairs and r — > ^v^Vt decays. The effect is largest (~ 5%) for cos 6 ^^ 
and large track momenta. 

Finally, hits are discarded from tracks where the corresponding hit from the following or the 
preceding track is missing. In those cases, the measured dE/dx is overestimated since the charge 
is not correctly distributed among the hits but is assigned to one hit only. By discarding this kind 
of hits, (3 — 5)% of all dE/dx hits are lost. 

The effect of all corrections can be seen in Figure |31 It shows the dE/dx pull distribution under 
a pion-hypothesis for all like-sign'' tracks from 3-prong tau decays. The solid line is obtained by a 
fit of three Gaussians to the measured pull distribution in 9 variable momentum bins from 3 GeV 
to 40 GeV. Since the number of muons in this sample is very small and their energy loss is very 
similar to the one for pions in the momentum range considered here, separate Gaussians for pions 
and muons are not needed. The free parameters in this fit are two of the three normalizations 
(the third one is constrained so that the sum of all three is equal to the observed number of 
tracks), one mean to allow for an overall shift (the position of the other two is calculated using 
the corresponding prediction from the Bethe-Bloch formula) and the width of the Gaussian which 
describes the contribution from pions. The deviation of this parameter from unity is used to 

^Tracks with the same charge as the initial r lepton. 
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Figure 3: Pull distribution obtained under a pion hypothesis for all tracks in 3-prong t lepton decays 
with a minimum momentum of 3 GeV and a minimum number of 20 hits in the dE/dx 
measurement. Only the tracks with the same charge as the decaying r lepton are shown. 
The solid points with error bars are data after all corrections and the function shows the 
expectation as explained in the text. The open points in the range between -8 to -2 show 
the same distribution hut without the corrections mentioned in the text. The smaller plot 
at the bottom shows the ratio of the full data points to the sum of the functions. 

obtain correction factors for the error of the energy loss measurement. The width of the two other 
Gaussians is calculated assuming that the relative error is constant. 

From the measured energy loss, its error and the expectation calculated using the Bethe-Bloch 
equation, x^ probabilities are calculated that the measured energy deposition is in accordance 
with the expectation for a given particle type. Pion- and kaon-weights, W-,^ and Wk, as used in 
this paper, are then calculated by taking one minus the value of this probability. These weights 
acquire a sign depending on whether the actual energy loss lies above or below the expectation for a 
certain particle hypothesis. This means that W-,^ is expected to be close to —1 for kaons since their 
energy loss per unit length is smaller in the momentum range relevant in this analysis. For electron 
tracks, W-,^ is expected to be close to +1 due to the higher energy loss in this case. Whenever these 
quantities are used in the selection, a cut on at least 20 dE/dx hits for this track is made implicitly. 



3.2 Photon Reconstruction and Identification of Neutral Pions 



The reconstruction of vr'^ mesons from photon candidates starts from an algorithm that has been 
used in previous OPAL publications (see e.g. [7]). It is based on the analysis of shower profiles in 
the ECAL as a function of the energy and direction of photons. In the fit each cluster is treated 



individually. Since the lead glass blocks in the barrel have a quasi-pointing^ geometry, only the 
lateral shower profile is considered. This algorithm provides the number of photon candidates and 
their corresponding 3- vectors that best describe the measured profile of the electromagnetic cluster. 
All photon candidates are assumed to originate from the interaction point (approximated by the 
primary vertex). Since the relative branching fractions of decays with and without tt^ mesons are 
different in strange and in non-strange final states, the parameters of the algorithm were optimized 
for this analysis as described below. 

Within each cluster only blocks with an energy of at least 150 MeV are considered and the total 
cluster energy must exceed 600 MeV. The expected mean energy deposited by a minimum ionizing 
particle is then subtracted from all ECAL blocks hit by a charged particle before the fit. If a photon 
candidate is too close to the entrance point of a track into the ECAL, the track's hadronic interaction 
can distort the photon energy measurement. Therefore, a minimum angle of 2.8° between a photon 
candidate and a track's ECAL entrance point is required. This value is obtained from studies of 



■'The longitudinal axis is pointing towards the interaction region. The blocks are tilted slightly from a perfectly 
pointing geometry to prevent neutral particles from escaping through the gaps between blocks. 
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Figure 4: Energy distribution of the reconstructed tt^ candidates. Plots (a) and (h) show the recon- 
structed energy for events with one track, plots (c) and (d) for events with three tracks. 
All plots are normalized to the number of t lepton decays. The arrows indicate the region 
selected. 
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the rate of fake tt*^ mesons in the decay r^ -^ K~Vr and subsequent optimization. To improve 
the energy resolution and the purity of the selection, a pairing algorithm is applied to recombine 
fake photon candidates with the closest photon candidates that were wrongly split up by the 
reconstruction procedure. The recombination is performed using the Jade jet finding scheme with 
the PO option [31]. Jet resolution parameter ycut values are optimized to obtain the best description 
of the number of expected photons using Monte Carlo events. Two ycut values are determined, one 
for clusters where a track is pointing to one of the blocks in the cluster and one for those without 
tracks. This is necessary since the hadronic interaction of charged particles disturbs the shower 
profile. The optimized ycut values are —3 and —4.6 for clusters with and without tracks, respectively. 
The angular resolution of this algorithm is 2.1° for clusters with tracks and 1.7° for clusters without 
tracks. Since these opening angles correspond directly to the energy of the neutral pion, photon 
candidates with an energy of more than 7.5 GeV are directly interpreted as neutral pion and the 
4- vector is corrected to account for the vr" mass. The energy of the reconstructed vr*^ candidates for 
events with one and three tracks is shown in Figure |1J The plots are normalized to the number of r 
decays in the event sample. Neutral pion candidates with an energy below 1.5 GeV in the 1-prong 
case and 2 GeV in the 3-prong case are rejected. 

For the remaining photon candidates, all two-photon combinations are tested. The combination 
which results in the maximum number of neutral pions with invariant two-photon masses not 
exceeding the vr'^ mass by more than 1.5(T is retained. A fit with a 7r° mass constraint is then applied 
to all vr'' candidates. All vr*^ candidates are assumed to originate from the primary interaction point. 



3.3 Identification of Kg 

The K*^ signal consists of 50% K^ and 50% Kg. The signature of a K^ decay is a large energy deposit 
in the hadron calorimeter without an associated track pointing to the cluster. The resolution of the 
OPAL hadron calorimeter would not allow for a clean reconstruction of this channel, thus it is not 
considered here. For the Kg, two decay modes are dominant, Kg -^ -k^-k^ (~ 32%) and Kg -^ 'k^tt~ 
{k, 68%). In this analysis only the latter is considered since the photon reconstruction algorithm 
exploits the quasi-pointing geometry of the electromagnetic calorimeter in the barrel (see Section 
?2|) . Thus only photons from the primary vertex can be properly reconstructed. 



The selection starts by combining each pair of oppositely charged tracks. Each track must have a 
transverse momentum with respect to the beam axis of px ^ 150 MeV, a minimum of 20 out of 159 
possible hits in CJ, at least 20% of all geometrically possible hits and a maximum x^ for the track 
fit of 50. For each combination of tracks, their intersection points in the plane perpendicular to 
the beam axis are calculated. The one with a radius less than 150 cm is selected as the secondary 
vertex. If two vertices are found that satisfy this condition, the one with the first measured hit 
closest to the intersection point is selected. In addition, the z-coordinate of the vertex has to satisfy 
|2;v| < 80 cm. 

The variables used in the selection described below are shown in Figure El For each Kg candidate, 
the angle between the reconstructed momentum of the candidate vertex and the Kg direction of 
flight must be less than |0v| < 0.5°. If the number of dS/dx-hits is more than 20 on at least one of 
the tracks, the energy loss information is used to veto photon conversions. All candidates where at 
least one track has a pion weight of more than 0.98 are rejected. To further reduce the background 
from photon conversions, the invariant mass of the Kg candidate is calculated using electron hy- 
pothesis for both tracks. All combinations with a mass rriee < 0.1 GeV are rejected. The remaining 
background that mainly comes from 3-prong r lepton decays or from wrong track combinations is 
rejected by applying cuts on the impact parameters (with respect to the nominal interaction point 
do) of the two tracks. The do values must have opposite sign according to the OPAL convention 
and the absolute values have to satisfy min(|(iQ |, I^q |) > 0.03cm and max(|(iQ |, |(iQ |) > 0.1cm. 
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Figure 5: Variables used in the Kg selection. A detailed description of all variables is given in the 
text. The dots represent the data and the open histogram, is Monte Carlo signal. The 
shaded areas show the background where photon conversions are marked separately. The 
arrows indicate the region selected. For all plots, all selection cuts have been applied 
except for the cut on the variable shown. All plots are normalized to the number of r 
decays. 
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Figure 6: Result of the Kg selection. Plot (a) shows the invariant mass distribution of the Kg 
candidates under pion hypothesis before the kinematic fit. Plot (b) shows the radius of 
the reconstructed secondary vertex and plot (c) the distribution of the x^ -probability of 
the 2C-fit. A cut is applied on the probability at 10"^. The dots represent the data and 
the open histogram is the Monte Carlo signal. The shaded areas show the background 
where photon conversions are marked separately. All plots are normalized to the number 
of T decays. 



The remaining Kg candidates must have a momentum of pj^o > 3 GeV. A 3D vertex fit is apphed to 

each candidate, that includes a constraint of the invariant two-track mass under the pion hypothesis 

to the nominal Kg mass. This is a 2C fit and a cut on the \^ probability at 10~^ is applied. 

The invariant two-track mass under the pion hypothesis before the kinematic fit can be found in 

Figure El together with the radius of the reconstructed secondary vertex and the x^ probability 

of the kinematic fit. In Figure |H{a) the vrvr-invariant mass spectrum is shown without any vertex 

constraint. If more than one Kg candidate shares the same track, the one with the smallest deviation 
„„ u„. 



from the nominal Kq mass before the fit is selected. 



After this selection procedure, a total of 535 Kg candidates remain with an estimated purity of 
82%. About 70% of the background consists of wrong combinations of tracks, and 30% comes from 
photon conversions. In one data event, two Kg candidates are found within one cone. This event 
is considered to be background. 
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4 Identification of r Final States 



For the selection of the various final states, a cut-based procedure is used where each r decay is 
treated independently. For all selected decay modes, the cone axis, calculated from the momenta 
of all tracks and neutral clusters identified in the electromagnetic calorimeter, must have a polar 
angle within | cos^| < 0.68 for the reasons explained above. Each selected cone must have at least 
one good track coming from the interaction point and the summed momenta of all tracks have to 
be less than the beam energy. Since there is at least one hadron in the final states considered here, 
the total energy deposited in the hadron calorimeter within the cone is required to exceed 1 GeV. 



4.1 Corrections to the Invariant Mass Spectra 

The spectral function from r lepton decays is a weighted invariant mass distribution of all hadronic 
final states with strangeness. The various final states have different experimental resolutions and 
migration effects that must be corrected for. To correct the observed data, the following method 
was used. The elements Cij of the inverse detector response matrix were determined directly from 
Monte Carlo. They represent the probability that an event reconstructed in bin j was generated 
in bin i. To calculate this probability Monte Carlo samples for the signal channels with mass 
distributions according to phase space were used. The corrected distribution was then obtained by 

Corrected _ V^ „ /^DATA _ Background n /r)\ 

J 

where g^ is the number of events in the data with a reconstructed mass in bin j, g- '^ s^°^'^ 

is the number of background events predicted in bin j and g9°^^'^'^^^'^ the number of events in mass 
bin i after correction. 

The corrected distributions will in general be biased towards the Monte Carlo input distributions. 
To reduce the bias from this approach, the method was applied iteratively. The result of the 
preceding iteration was used to refine the elements of the inverse detector matrix. The optimal 
number of iterations was determined, using Monte Carlo simulations, to be two for all final states 
considered here. The corresponding systematic uncertainties are discussed in Section 15.41 Finally, 
detection efficiency corrections were applied to the corrected invariant mass distributions. The 
bin width of 150 MeV of the measured invariant mass spectra were chosen according to the mass 
resolution of the r decay channels measured. 



4.2 (K^) Final States 

The (Ktt)" mass spectrum consists of two measured modes, ¥^~'k^Vt and ¥>PiT~Vr- From the latter 
decay mode, only the decays Kg -^ vr+vr" are measured. The K*(892) dominance is well established 
in both cases. 



4.2.1 K-TT^Vr 

In the K~-K^iyr selection, exactly one good track coming from the primary vertex is required. This 
track must have a minimum momentum of p > 3 GeV. For the track to be selected as a kaon, the 
pion weight has to satisfy Wt^ < —0.98 and the kaon weight Wk < 0.6. Furthermore, exactly one 
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identified neutral pion is required with E^^^ > 1.5 GeV selected using the algorithm explained in 
Section [3.21 The variables used in this selection are shown in Figures |S1 and [71 



From this selection 360 events are seen in the data with 190.5 background events predicted from 
Monte Carlo. This corresponds to a background fraction of 54%. The invariant mass spectrum 
can be seen in Figure [HI The mass resolution in this channel is approximately 40MeV. The 
main background comes from t~ — > K~i/^ decays where one fake neutral pion was reconstructed. 



Additional sources of background are r 



vr 



'ir^Ur, 



T ^ K K'^z^t- and r — > K^K tt'^z^t- where 



the K'^ is a Kg decaying to two neutral pions or a K^ which does not decay within the jet chamber. 
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Figure 7: The pion weight W-,^ for 1-prong t decays in the range (—1 : —0.8). The dots are the data 
points. The open (shaded) histogram denotes the contribution from kaon (pion) tracks 
as predicted by the Monte Carlo simulation. The plot is normalized to the number of r 
decays. Events on the side of the direction of the arrow are considered kaon candidates. 
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Figure 9: Result of the K~7r'^zvT- selection. Plot (a) shows the measured invariant mass spectrum. 
The dots are the data, the open histogram is the Monte Carlo signal and the shaded area 
is the background. Plot (b) shows the background subtracted spectrum, (c) shows the 
corrected spectrum. Note that the correlation between neighboring bins in the corrected 
spectrum is of the order of 50%. Plot (d) is the selection efficiency as a function of the 
invariant mass. 



4.2.2 KO 



vr v-r 



The selection is very similar to that for the KTit^Vt final state. Here exactly one identified Kg is 
required using the procedure from Section [3.31 In addition, one good track from the primary vertex 
is required. If the momentum of this track lies above the kinematically allowed minimum for a kaon, 
the same identification procedure as mentioned above is applied to veto decays r~ — > K~K'^i/,-. 
Only events with zero reconstructed vr*^ mesons are accepted. The variables used in the selection 
are shown in Figure [TUl 

From this selection 361 events are expected with a background fraction of 47%, and 344 are seen in 
the data. The main background contributions come from decays t~ -^ K'^K T^'fT, t~ — > K~K'^i/t- 
and T~ — > }^^it~'k^Vt- where the neutral pion escapes detection. The invariant mass spectrum for 
this channel is shown in Figure ITTl The mass resolution in this channel is approximately 60MeV. 
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4.3 (Ktttt) Final States 

The (Ktttt)^ final state consists of the decay modes: K~ n^ tt^ Ur , K^ir^ir^h'r and K^tt^tt^i't-. To 
select these final states, the following procedure is applied. 



4.3.1 K^TT+vr^i/^ 

The selection starts by requiring exactly three good tracks coming from the interaction point. These 
tracks are fitted to a common vertex and the fit probability is required to be larger than 10^'''. In 
addition, each pair of oppositely charged tracks has to fail the selection criteria for neutral kaons 
as defined in Section f3. 31 These two requirements reduce the background from photon conversions 
and decays containing Kg. 

To identify the kaon in the decay, one of the like-sign candidate tracks must have p > 3 GeV and 
Wtj- < —0.9. To further reduce the pion background among these candidate tracks, Wk and Wt^ 
are input to a neural network. The track is rejected if the output of the neural network is below 
0.3 (see Figure IT^ c)). Exactly one like-sign track is allowed to fulfill these requirements, otherwise 
the decay is treated as background. If the momentum of the unlike-sign track is consistent with 
the T~ — > TT'K'K^Vr hypothesis, this r decay is only accepted if W-n- > —0.95 (see Figure El^d)). 

The algorithm for identifying neutral pions (Section 13. 2|) is then applied to the selected cones. For 
this channel, the number of reconstructed vr" mesons with an energy greater than 2 GeV is required 
to be zero (see Figure IT^ b)). Otherwise this r decay is treated as background. To further improve 
the purity of the selection, the cosine of the decay angle in the rest frame of the r lepton, the so- 
called Gottfried- Jackson angle Q* is calculated. The cos0* distribution is shown in Figure IT^ al. 
For events where the kaon hypothesis was applied to the wrong track or the number of identified 
neutral pions does not correspond to the true number, this calculation leads to unphysical values of 
that variable. Due to resolution effects, correctly identified signal events can also give values beyond 
±1. Therefore, a cut was applied at cos 0* = ±1.2. The contribution from t~ — > K^7r+7r~(n7r'^)i/,- 
events is included in the background estimate. 

From this selection 269 events are seen in the data with a contribution of 149.8 background 
events predicted from Monte Carlo. This corresponds to a background fraction of 63%. The 
main background contribution comes from decays r~ — > tt~ tt^ tt^ Vr , t~ — > 7r~K~K+i/,- and 
T~ -^ K~ TT^ TT~ TT^ Vr , where the vr'^ meson escapes detection. The invariant mass spectrum can 
be found in Figure IT^ The mass resolution in this channel is approximately 20MeV. 



4.3.2 K^TT-ir^Ur 

Exactly one identified Kg and exactly one tt^ meson with an energy E'^ > 2 GeV is required for 
this final state. The pion candidate track has to satisfy the same requirements as for (Kvr)" final 
states. 

From this selection, 65 events are expected and 67 are seen in the data with a background fraction 
of 72%. The main background contribution comes from decays r~ -^ K^-K^Ur where the neutral 
pion escapes detection. The invariant mass spectrum can be seen in Figure ITU The mass resolution 
in this channel is approximately 100 MeV. 
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Figure 12: Variables used in the K~7r+7r~z/T- selection. Plot (a) shows the cosine of the Gottfried- 
Jackson angle Q* , (b) the number of reconstructed neutral pions and (c) the output of the 
neural network. The dots are the data, the open histogram is the background prediction 
from the Monte Carlo and the shaded histogram denotes the contribution from K^vr+vr^n 
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Figure 13: Result of the K^tt^tt^Ut selection. Plot (a) shows the measured invariant mass spec- 
trum. The dots are the data, the open histogram is the Monte Carlo signal and the 
shaded area is the background. Plot (h) shows the background subtracted spectrum, plot 
(c) is the corrected spectrum. Note that the correlation between neighboring bins in 
the corrected spectrum is of the order of 50%. Plot (d) is the selection efficiency as a 
function of the invariant mass. 
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Figure 14: Result of the K^tt^tt^Vt selection. Plot (a) shows the measured invariant mass spectrum. 
The dots are the data, the open histogram is the Monte Carlo signal and the shaded area 
is the background. Plot (b) shows the background subtracted spectrum, (c) the corrected 
spectrum. Note that the correlation between neighboring bins in the corrected spectrum 
is of the order of 50%. Plot (d) is the selection efficiency as a function of the invariant 
mass. 
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4.4 (KyrTTTr) Final States 



The (Kvrvrvr)" signal consists of the following final states: K~ n'^ tt^ n^ i^r , K^vr^vr^vr^i/. 

and K^tt^tt^tt^Ut 

investigated. 



0^0,, K-TrOvrOvrOi/, 
From these, only the first one which has the highest branching fraction is 
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Figure 15: Result of the }i~ tt^ tt^ t:^ u-r selection. Plot (a) shows the measured invariant mass spec- 
trum. The dots are the data, the open histogram is the Monte Carlo signal and the 
shaded area is the background. Plot (h) shows the background subtracted spectrum. 

The same procedure as for the K^tt+tt^z^t- channel is used. In addition, one identified vr*^ meson 
with an energy of more than 2 GeV is required. The invariant mass spectrum can be seen in Figure 
1151 The mass resolution in this channel is approximately 60MeV. From this selection, 14 events 
are seen in the data with a contribution of 10 events from background. The selection efficiency is of 
the order of 1%. The main background contribution comes from r~ — > K~7r^7r^^'T- decays, where 
one fake neutral pion was identified. 

Since the number of signal events in this final state is not significantly different from zero, this 
channel is not considered any further in this analysis. For the spectral function, the Monte Carlo 
prediction has been used instead. 



5 Results 



5.1 Branching Fractions 



The measured data used in the spectral function analysis allows the determination of competitive 
branching fractions for the channels r^ -^ Y^^ir^Vr and r^ -^ K^ir^ir^Ur. The branching fractions 
are determined in a simultaneous x^ fit? taking all measured final states into account. The fit 
function is 



AT, = N^on-r + (1 _ f 



non- 
bkg 



N'2^e.,BjF^ 



Bias 



where i is the signal channel under consideration and index j runs over all r decay channels. Ni 
is the number of decays observed in the data, i.e. the number of r cones passing the cuts. The 
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number of decays from the non-r background is denoted by N^°^ ^, /^^g ^ is the fraction of non- 
r background, N"^ is the number of r events, £ij the efficiency matrix, F^^^ the bias factor for 
the respective channel due to r preselection cuts and Bj are the fitted branching fractions. The 
branching fractions for the two signal channels were allowed to vary freely, the branching fractions 
of all other r decay channels contributing to the background prediction were allowed to vary only 
within their PDG errors [14]. The selection quantities used in the fit together with the contributions 
from the background channels and their branching fractions used in the Monte Carlo simulation 
are shown in Table |21 From the fit, the following results are obtained: 

B{t- -^ K"7r°i/^) = (0.471 ± 0.064stat ± 0.022sys)% 
B{t- ^K-Tr+ IT- {OTT^,ex.K^)iyr) = (0.415 ± 0.059stat ± 0.031sys)%, 

where the possibility that the vr+vr" pair in the K~ n'^ n" i^r final state comes from a Kg has been 
excluded. For the estimation of the systematic uncertainty, the following sources are considered. 
They are summarized in Table El The total systematic uncertainty is obtained by adding the 
individual contributions in quadrature. 

• Energy loss measurement (^dE/dx)' 

In the selection, the specific energy loss dE/dx is used to separate pions from kaons. Cuts on 
corresponding weights are applied which are calculated from the pull distribution (see Figure 
EJ. A possible shift in this quantity can lead to a systematic misidentification of tracks. The 
pull distribution is therefore shifted within the error on its mean and the selection procedure 
is repeated. The difference between the branching fractions obtained with and without the 
shift applied is the systematic uncertainty. 

• Energy scale in vr*^ reconstruction (A^;): 

The energy resolution can be tested by measuring the invariant two-photon mass from vr*^ 
decays. A systematic shift in the observed mass in the data compared to the detector simu- 
lation can be translated into a scale factor for the reconstructed photon energies. Deviations 
of (0.5 lb 0.9) MeV from the nominal vr*^ mass have been observed [7], corresponding to a 
scale factor of 1.004 it 0.007. The energies of the reconstructed photons in the Monte Carlo 
samples were therefore varied by ±0.7%. The difference between the branching fractions 
obtained with and without the variation applied is the systematic uncertainty. 

• Momentum scale (Ap): 

The systematic uncertainty connected with the momentum scale was tested using Z^ — > fJ-~ n'^ 
events [7]. The difference in momentum resolution between data and Monte Carlo as a 
function of cos 6 was studied. To assess the systematic uncertainty Ap in hadronic r decays, 
all particle momenta in the Monte Carlo were varied accordingly. The difference in the result 
with and without this variation is quoted as a systematic uncertainty. 

• Monte Carlo statistics (Amc): 

The precision of the background prediction depends on the Monte Carlo statistics used in the 
selection procedure. Therefore, the number of background events selected is varied randomly 
within its statistical uncertainty. The observed spread in the branching fraction due to this 
variation is quoted as a systematic uncertainty. 

• Bias factor (A^Bias): 

The bias factors determined from the Monte Carlo are varied by their uncertainty and the 
branching fractions are then refitted. The observed spread due to this variation contributes 
to the total systematic uncertainty. 
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T- -^ K-TtO//^ 








No. of Events 


360 




Selection Efficiency /% 


8.42 ± 0.17 








Preselection Bias Factor 


1.016 ± 0.011 








Non-r Background Fraction 


0.006 ± 0.004 








r Background Fraction 


0.540 ± 0.027 








13.5% 


0.051 ± 


0.005 


25.41 ± 0.14 


K-KOyrOj^^ 


9.9% 


6.0 ± 


0.3 


0.155 ± 0.020 


K~Vt- 


8.1% 


1.25 ± 


0.06 


0.686 ± 0.023 


K-ROzv^ 


7.3% 


4.5 ± 


0.2 


0.154 ± 0.016 


Tr^vrOvrOz., 


6.2% 


0.07 ± 


0.01 


9.17 ± 0.14 


K-ttOttOz., 


5.0% 


9.6 ± 


0.5 


0.058 ± 0.023 


K-TrOTrOvrOj., 


2.6% 


8.9 ± 


0.6 


0.037 ± 0.021 


other 


1.4% 










Bkg. Fraction 


Efficiency /% 


^PDG/% 



r ^ K vr+vr i 


^T 






No. of Events 


269 




Selection Efficiency/% 


6.59 ± 0.06 






Preselection Bias Factor 


0.953 ± 0.013 






Non-r Background Fraction 


0.007 ± 0.006 






r Background Fraction 


0.631 ± 0.044 






21.6% 


0.15 ± 0.02 


9.22 ± 0.10 


TT-K-K+Ur 


10.3% 


3.9 ± 0.2 


0.161 ± 0.019 


TT^TT^TT^TT^Ur 


8.1% 


0.5 ± 0.1 


4.24 ± 0.10 


K^TT^ir^-K^VT- 


6.7% 


2.7 ± 0.2 


0.064 ± 0.024 


other 


16.4% 








Bkg. Fraction 


Efficiency /% 


^PDG/% 



Table 2: Quantities used in the fit for the branching fractions. The errors quoted for efficiency, bias 
factor and background fractions are from Monte Carlo statistics only. The last column 
contains the branching fractions for the background channels used in the Monte Carlo 
simulation [I4]. 





T~ -^ K-TT^Ur 


r ^ K TT~^TT V-T 


Energy Loss Measurement (dE/dx) 
Energy Scale 
Momentum Scale 
MC Statistics 
Bias Factor (F^^^^) 


0.012 
0.010 
0.003 
0.014 
0.004 


0.019 
0.011 
0.003 
0.021 
0.005 


Total 


0.022 


0.031 



Table 3: Individual contributions to the systematic uncertainty of the branching fraction measure- 
ments as explained in the text. 
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5.2 Improved Averages for B{t 



K 



-TrV. 



and B(t 



K- 



TT^TT U^ 



For the determination of the spectral function and its moments described below, new average values 
for the branching fractions of the decays r~ -^ K'tt^u-j- and r~ — > K^ tt^ -k^ u-j- are determined. In 
addition to the results obtained here, the same measurements are used as inputs for the calculation 
as in [14]. For the channel t~ — > K~7r^7r~^'T-i the previous result from OPAL [32] was replaced by 
the measurement obtained in this analysis. Also, the CLEO result was updated using [33]. The 
new averages are: 



(0.453 ±0.030)% 
(0.330 ±0.028)%. 



5av(T- ^ K-vrV,) 

The measurements used, together with the averages given in [14] and the improved values for the 
branching fractions obtained here are displayed in Figure [TBI 

The branching fraction obtained from this analysis for the K'^tt^i't- channel is consistent with the 
previous measurements within the errors quoted. The value obtained for the K~7r^7r~z^T- channel 
is consistent with the new measurement from CLEO and the theoretical prediction in [22]. The 
ALEPH [34] result differs from these values by roughly 2.5a. 
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Figure 16: New averages for the branching fractions for r^ -^ K^tt^i't- (a) and t~ -^ K^vr^vr^i/,- 
(h). The previous world average as given by the PDG is indicated by the shaded band. 
The measurements used in the calculation of the weighted mean and the new averages 
are shown. OPAL03 denotes the values for the branching fractions obtained from this 
analysis. 



5.3 The Strangeness Spectral Function 



The hadronic decay of the r lepton is commonly written in terms of the so-called spectral functions 
ui(s), ao/i(s) for the non-strange part and ■Vg/iC'^) ^'^d 0^,-^(5) for the strange part. The functions 
V and a are the vector (V) and the axialvector (A) contributions, respectively, while the subscript 
denotes the angular momentum J . The variable s is the invariant mass squared of the hadronic 
system. The spectral function is experimentally determined by measuring the invariant mass spec- 
tra of the given hadronic modes and normalizing them to their respective branching fractions. The 
contributions to the total strangeness spectral function then read: 
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'VeVr) 
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V/A 



ds 
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B{T^{V/A)^^=-^^J=^)vr) 1 diVv/A 
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e VpVt 



N^ 



V/A 



ds 



(3) 



(4) 
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where |Kis| = 0.2196±0.0023 [14] is the CKM weak mixing matrix element, m^ = (1 776.9lo;27) MeV 
[35] and Sew = 1.0194 it 0.0040 [36] is an electroweak correction factor. The total strangeness 
spectral function (v + a) is then obtained by adding the individual contributions. To disentangle 
the vector and the axialvector contributions for the spin-1 part, a detailed analysis of the resonance 
structure of the measured spectra is necessary which is not done here due to the limited statistics. 
The kaon pole contributes to the pseudoscalar spin-0 part. 

Summing the individual contributions (see Equations El and ^ , the strangeness spectral function 
is obtained. The Monte Carlo prediction of the total strangeness spectral function as a function of 
the invariant mass squared is displayed in Figure 1171 The improved version of the r Monte Carlo 
as explained in Section [2.31 has been used here. The branching fractions as given in Tabled have 
been used to weight the individual mass spectra. For illustration purposes the spectral function is 
shown using two different binnings. A non-equidistant binning is chosen which corresponds to a 
bin width of 50 MeV and 150 MeV in the invariant mass, respectively. The errors given here are 
from Monte Carlo statistics only. 

The spectral function obtained from the data is displayed in Figure ITHl For t~ -^ K^tt'^z^t- and 
r^ -^ K^tt^tt^^'t-, the new average branching fractions and their respective errors as given in 
Section [5.21 are used. The binning chosen in this plot corresponds to a bin width of 150 MeV in the 
invariant mass and is governed by the mass resolution of the K^vr"?!'^!/,- final state. The dots with 
error bars represent the inclusive spectrum. The inner error bars are the statistical uncertainties. 
They include the uncertainty on the efficiency and on Monte Carlo statistics. The total error is 
calculated by adding up the statistical and systematic uncertainties (explained in the next section) 
in quadrature. The numerical values are given in Table |3 The systematic uncertainty is dominated 
by the uncertainty on the r branching fractions. 

For the (Kvr)" final state both r decay channels K^tt~i't- and K~7r^i/^ are measured. The channel 
K^ijUr which also contributes to the two meson final state is taken from Monte Carlo. For the 
(Ktttt)^ final state the spectra K.^tt'tt^i't and K~7r"'"7r~i/T- are measured. The contribution from 
the decay K'tt^tt^Ut- is added from Monte Carlo as well as the K'rjir^Ur channel which also con- 
tributes to the three meson final states. For the (KTrvrvr)" spectrum, which consists of the channels 
K~7r"'"7r~7r''^'T-, K^tt~tt^tt^Ut-, K~Tr^'ir^TT^h'T- and K'^Tr^vr+vr"^',-, the prediction from the Monte Carlo 
is taken. 



5.4 Systematic Uncertainties on the Spectral Function 

The sources for possible systematic uncertainties listed below have been considered. Since the 
individual contributions are different for the different final states, the error is given for each bin in 
s separately in Table 01 

• Energy loss measurement {A.^E/dx)' 

The systematic variation is described in Section [5. II In addition, the momentum dependence 
of this variation has been investigated. No significant infiuence has been found. 

• Energy scale in vr*^ reconstruction (Ae): 

The systematic variation is described in Section [5.11 

• Momentum scale (Ap): 

The systematic variation is described in Section [5.11 

• PDG errors on branching fractions (Ab): 

The dominant contribution to the systematic uncertainty comes from the uncertainty in the 
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r branching fractions (see Table ^ . The branching fractions are varied randomly and the 
difference in the result is quoted as systematic uncertainty. The channels which populate the 
region of high s have branching fractions with relative errors close to 100% leading to large 
uncertainty in the spectral function itself. This also covers the uncertainties on the shape. 

Kg identification (Aj^o): 

A possible origin for systematic effects in the Kg identification is the estimation of the back- 
ground using Monte Carlo. In particular the number of photon conversions found in r lepton 
decays is not perfectly modeled. The cut on the x^ probability of the 2C constrained fit to Kg 
(see Section [3. 3|) has been varied from 10^^ to 0.01 to estimate a possible systematic effect. 
This cut reduces the number photon conversions in the sample by one order of magnitude. 

Mass correction procedure (Amcorr): 

In order to assess the systematic uncertainty associated with the mass correction procedure, 
two possible sources have to be considered. The effect due to different input mass spectra 
and the effect due to the choice of the number of iterations. They have been studied using 
high statistics event samples. The systematic deviations using either flat, phase space, or 
resonance shaped input spectra (with an optimized number of iterations for each scenario) 
are about 5% in each mass bin. The deviations are largest when a flat input distribution 
is assumed. Because the true mass spectrum is certainly not flat, we consider 5% to be a 
conservative estimate for Amcorr- The effect of the number of iterations was estimated using a 
phase space distribution for input spectrum and performing one additional iteration step. The 
change was found to be negligible compared to the effect of varying the input distribution. 



(s-range)/GeV2 


Ab 


^dE/dx 


AKO3 


Ae 


Ap 


^mcorr 


A tot 

sys 


Astat 


V + a 


(0.53,0.77) 


0.04 


0.006 


0.006 


0.007 


0.003 


0.06 


0.07 


0.17 


1.17± 0.18 


(0.77, 1.06) 


0.13 


0.011 


0.011 


0.014 


0.001 


0.11 


0.17 


0.18 


2.27± 0.25 


(1.06,1.39) 


0.08 


0.003 


0.003 


0.004 


0.001 


0.03 


0.09 


0.07 


0.69± 0.11 


(1.39,1.77) 


0.18 


0.005 


0.005 


0.005 


0.002 


0.05 


0.18 


0.19 


0.90± 0.26 


(1.77,2.19) 


0.32 


0.006 


0.007 


0.007 


0.003 


0.06 


0.33 


0.25 


1.22± 0.41 


(2.19,2.66) 


0.35 


0.007 


0.009 


0.009 


0.003 


0.07 


0.36 


0.49 


1.44± 0.61 


(2.66,3.17) 


0.30 


0.007 


0.008 


0.008 


0.003 


0.07 


0.31 


0.85 


1.35± 0.90 



Table 4: Result for the strangeness spectral function. The table shows the values of the total 
strangeness spectral function {v + a) and the statistical and the total systematic uncer- 
tainties for every bin in s. The individual contributions to the systematic uncertainty are 
discussed in the text. The total uncertainty quoted is the quadratic sum of the statistical 
and systematic uncertainties. 
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5.5 The Spectral Moments k;[ 



The moments of the spectral function are defined as: 



R 



kl 

T,S 



m^ 



d. 1-^ ^ 



S(r^(y/A)(^=-i'-^=o/i)i/^) 1 diVy/A 



(5) 
/o \ '"-t/ V'"-?/ B{t ^ e-UeVr) ^V/A ds 

where the sum runs over all final states with open strangeness. The values measured for the mo- 
ments kl = (00, 10, 11, 12, 13, 20, 21, 30, 40), their statistical and systematic uncertainties are 
given in Table El The various sources contributing to the systematic uncertainty have been dis- 
cussed in Section 15.41 The value for T&^g is calculated from the branching fractions only and is 
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therefore independent of the measured spectra. In addition, the CKM-weighted difference of the 
corresponding strange and non-strange moments is given: 



5R 



kl 



R: 



kl 



iKdp 



nkl 

\V |2 
''lis 



(6) 



(see Table ISJ where R^'-g are the strange moments and R^^non-s ^^ *^^ ^^™ '-'^ ^^^ vector- and 
axialvector non-strange moments from [7], using updated branching fractions. The values of the 
matrix elements used for the non-strange and strange moments are iKidl = (0.9734 it 0.0008) and 
iKisI = (0.2196 lb 0.0023), respectively [14]. Systematic errors which are common to both, strange 
and non-strange moments, such as the energy scale error and the momentum scale for tracks, are 
100% correlated and are treated accordingly. 
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Figure 17: The strangeness spectral function as determined from Monte Carlo. For illustration 
purposes two different binnings were used - (a) 150 MeV and (h) 50MeV in mass. The 
dots represent the inclusive spectrum. The white histogram denotes the contribution from 
the (K7r)~ and (Kr/)~ final states, the dark shaded histogram those from the (KvrTr)^ and 
(K?77r)^ channels and the light shaded area shows the contribution from (Kvrvrvr)^ final 
states. The leftmost bin represents the kaon pole which is the only contribution to the 
spin-0 part. The errors shown include the statistical uncertainty of the Monte Carlo as 
well as the uncertainties on the branching fractions. 
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Table 5: The spectral moments for kl = (00,10,11,12,13,20,21,30,40). The tables include the values for the non-strange moments, the strange 
moments, the weighted difference and the quotient of strange to non-strange moments. The values for the non-strange moments given are 
calculated from [7] using updated branching fractions. For the strange moments, statistical and systematic uncertainties are listed separately 
for the individual sources. The statistical uncertainty also contains the uncertainty on the branching fractions. The uncertainty on R^^ is 
the total error from the branching fractions as given in [14]. On the right hand side of each table, the correlations are given in percent. 
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Figure 18: The spectral function from strange r lepton decays. The dots show the inclusive spectrum 
as measured in OPAL. The white histogram shows the exclusive (Kvr + K?])" spectrum, 
the dark shaded area the (Kvrvr + Kr/vr)" contribution and the light shaded area the con- 
tribution from the (KvrTrvr)^ final states. The error bars show the statistical uncertainty 
(inner error bar) and the systematic uncertainty added in quadrature (total error). 
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6 Conclusions 



The spectral function in r lepton decays into vector and axial vector final states with open 
strangeness, as well as the spectral moments R^^g, have been measured with the OPAL detec- 
tor, operating on the Z^ resonance. The selection of the decay channels requires a detailed 
study of the charged kaon identification using dE/dx, neutral pion reconstruction in the elec- 
tromagnetic calorimeter, and identification of Kg. For the decays t~— > (Kvr)"!/,-, (Kvrvr)"!/,- and 
(KvrTrvr)"!^,-, 94% of all final states have been experimentally accessed; for the remaining 6.6%, 
as well as for the strange final states including r] mesons, Monte Carlo simulations have been 
used. The reconstructed strange final states corrected for resolution effects and detection effi- 
ciencies yield the strangeness spectral function of the r lepton. The spectral moments -R^ 5 for 
kl = {10, 11, 12, 13, 20, 21, 30, 40}, as well as the ratio of strange to non-strange moments, up- 
dated from [7], and the weighted differences SR!^ are given, which are useful quantities for further 



theoretical analysis. For the decay channels r 
branching fractions are obtained: 



K tt"iJt- and r 



K TT'^TT Uj- competitive 



B(t- 



K 



B{t- 



-TT^. 



K vr^vr u-r- 



(0.471 ± 0.064stat ± 0.022 
(0.415 ± 0.059stat ± 0.031 



sys; 



sys; 
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